August, 1999 
OU-HEP-324 



Do we expect light flavor sea-quark asymmetry also for 
the spin-dependent distribution functions of the nucleon ? 

M. WakamatsuQ 

Department of Physics, Faculty of Science, 
Osaka University, 
Toyonaka, Osaka 560, JAPAN 

T. WatabeQ 

Research Center for Nuclear Physics (RCNP), 
Osaka University, 
Ibaraki, Osaka 567, JAPAN 



PACS numbers : 12.39.Fe, 12.39.Ki, 12.38.Lg, 13.40.Em 



Abstract 



After taking account of the scale dependence by means of the standard DGLAP evolution 
equation, the theoretical predictions of the chiral quark soliton model for the unpolarized 
and longitudinally polarized structure functions of the nucleon are compared with the recent 
high energy data. The theory is shown to explain all the qualitative features of the exper- 
iments, including the NMC data for -Ffl^) ~ -^r(^)) ^2 i^) / -^2 i^) ^ the Hermes and NuSea 
data for d{x) — u{x), the EMC and SMC data for g\{x), gi{x) and gfix). Among others, 
flavor asymmetry of the longitudinally polarized sea-quark distributions is a remarkable pre- 
diction of this model, i.e., it predicts that Ad{x) — Au{x) = C x°'[d{x) — u{x)\ with a sizable 
negative coefficient C ~ —2.0 (and a ~ 0.12) in qualitative consistency with the recent 
semi-phenomenological analysis by Morii and Yamanishi. 



^Email : wakamatu@milio.rcnp.osaka-u.ac.jp 
^Email : watabe@rcnp.osaka-u.ac.jp 



Undoubtedly, the most natural and widely-accepted explanation of the light flavor sea- 
quark asymmetry in the nucleon, revealed by the NMC measurements [|l|, would be the one 
by the pion cloud effects Accepting its essential validity, the NMC observation may be 

taken as the flrst clear evidence manifesting the nonperturbative QCD dynamics of spontaneous 
chiral symmetry breaking in high-energy deep-inelastic observables. This flavor asjTiimetry 
of the spin-averaged sea-quark distributions could have been predicted by some theoretical 
models which properly incorporate the chiral symmetry, but unfortunately, most theoretical 
explanations were given only after the NMC observation. 

Now, a natural question is "Do we expect light flavor sea-quark asymmetry not only for 
the unpolarized distributions but also for the spin-dependent ones?" In our opinion, to make a 
trustworthy prediction for spin-dependent sea-quark distributions, the chiral symmetry would 
be a minimum ingredient to be incorporated into model construction. Otherwise, it would 
be hard to explain the already-established flavor asymmetry of the unpolarized sea-quark 
distributions. Also desirable is a reasonable reproduction of the existing high energy data for 
the unpolarized and the longitudinally polarized structure functions. 

We claim that the chiral quark soliton model (CQSM) of the nucleon pf- \V2\ fulfllls the 



above requirements and that it is in this sense qualifled as a good candidate which is able 
to give some deflnite prediction for the flavor asymmetry of the spin- dependent sea-quark 
distributions. In the present report, we shall flrst compare the theoretical predictions of the 
CQSM with the recent high energy data for the unpolarized as well as the longitudinally 
polarized structure functions, including the NMC data for -Ffl^) ~ -^2^(2^)) F2{x)/F2{x) [|I|, 



the HERMES ll and NuSea (E866) data |T1 for d{x)-u{x), the EMC |Ti,[lT3,lT§ and SMC 
data for gfix), gi{x) and gf{x). After conflrming that the theory is able to explain all the 
qualitative feature of these experimental data, we flnally give a prediction for the longitudinally 
polarized anti-quark distributions Ad{x) — Au{x) in the proton. 

The chiral quark soliton model (CQSM) is a very simple model of the nucleon which 
maximally takes account of the spontaneous chiral symmetry breaking of the QCD vacuum. 
It is specifled by the effective lagrangian : 

£o = ^(i^-Me*^^'^-^(^)/-^- )?/^, (1) 

which describes the effective quark flelds with a dynamically generated mass M, interacting 
with massless (or nearly massless) pions. The theory is not a renormalizable one and it is 
deflned with some ultraviolet cutoff A. The above effective lagrangian was originally introduced 
by Diakonov and Petrov on the basis of the instanton-liquid picture of the QCD vacuum 
pO| . In this derivation, the dynamical quark mass M and the physical cutoff A can both be 
estimated as functions of the average instanton size p with 1/p — 600 MeV and the average 
separation between instantons R with 1/^ ~ 200 MeV. Their estimate gives M ~ 350 MeV 
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and A ^ 1/p ~ 600 MeV. To be more precise, the predicted dynamical quark mass M is 
momentum dependent, but in its application to the chiral soliton problem, it is customary to 
approximate it by a momentum independent constant. At the same time, the ultraviolet cutoff 
is fixed by the condition that the effective meson lagrangian derived from (|1]) reproduces the 
correct normalization of the pion kinetic term. For instance, in the Pauli-Villars regularization 
scheme, which is used throughout the present analysis, what plays the role of a ultraviolet 
cutoff is the Pauli-Villars mass Mpy obeying the relation : 



47r2 V M 

with the pion weak decay constant. Using the value of M ~ 375 MeV, which is favored from 
the phenomenology of nucleon low energy observables, this gives Mpy — 562 MeV in qualita- 
tive consistency with the estimate based on the instanton picture [^. (We refer to [|I^ for the 
justification of the single-subtraction Pauli-Villars scheme in the studies of quark distribution 
functions.) Since we are to use these values of M and Mpy, there is no free parameter addi- 
tionally introduced in the calculation of distribution functions. The model is also compatible 
with the idea of large QCD advocated by Witten many years ago |2^. According to him, 
in the limit of Nc = oo, the QCD is equivalent to an effective theory of mesons, and a baryon 
appears as a topological soliton of this effective meson lagrangian. Though an appealing idea, 
a practical problem of this scenario is that the relevant effective meson theory may not be 
so simple as that of the Skyrme model. In fact, who can imagine an effective meson theory, 
which is able to describe deep-inelastic scattering observables of the nucleon? On the contrary, 
the chiral quark soliton model is an effective theory, which realizes the idea of large Nc QCD 
in more economical way. In fact, in the stationary-phase approximation, which is exact in 
the large Nc limit the theory is known to have a self-consistent soliton-like solution of 
hedgehog shape. The energy degeneracy of this soliton configuration under the spatial and 
isospin rotation induces spontaneous rotational motion of hedgehog soliton. The semiclassi- 
cal quantization of this collective rotation leads to a systematic method of calculation of any 
nucleon observables, which is given in a perturbative series in the collective angular velocity 
operator Q. (This takes the form of a 1/Nc expansion, since Q itself is an 1/Nc quantity.) The 
numerical method of calculation of the static nucleon observables was established many years 
ago 0], and it has been successfully applied to many interesting observables However, 
the calculation of quark distribution functions, especially with full inclusion of the vacuum 
polarization contributions, is quite involved. The calculation of all the combinations of the 
twist-2 distributions has been completed only recently [Q, |TD[-|T2|. 

We first summarize in Fig.l our final results for the unpolarized and longitudinally polarized 
distribution functions of the leading twist 2. (In the present investigation, some refinement 
has been achieved in the numerical method of calculation of these distribution functions in the 
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small X region. As a consequence, there are some minor changes from the results reported in 
our previous paper [0, especially for the distribution Au(a;) + ^d{x). However, these changes 
are so small that they hardly leave appreciable influence on the structure functions obtained 
after Q^-evolution.) 



Figure 1: The theoretical predictions of the CQSM for the unpolarized distributions u{x)+d{x) 
and u{x) — d{x) as well as for the longitudinally polarized distributions Au{x) + Ad{x) and 
Au{x) — Ad{x). In all the figures, the long-dashed and dash-dotted curves respectively stand 
for the contributions of the discrete valence level and that of the Dirac continuum in the 
self-consistent hedgehog background, whereas their sums are shown by the solid curves. The 
distribution functions in the negative x region are to be interpreted as the antiquark ones 
according to the rules (0)^(0) • 

We recall that, within the theoretical framework of the flavor SU(2) CQSM, any nucleon 
observables are evaluated in either of the isoscalar combination or the isovector one In 
fact, theoretical formulas for them have totally dissimilar form reflecting the fact that they 
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have different dependence on the collective angular velocity f2 given as p,[|Tl|],||T2| 

u{x) + d{x) ~ A^,0(^]°) ~ 0{N^), (3) 

u{x)-d{x) ~ NcO{Q^) ~ 0{N^), (4) 

AM(x) + Arf(x) ~ iV,0(l]^) ~ 0(Ar°), (5) 

Au{x)-Ad{x) - iV,[0(fi°) + 0(f]i)] - 0(iV]) + 0(A^°), (6) 

where use has been made of the fact that Q scales as l/Nc- In Fig.l, quark distributions with 
negative x are to be interpreted as antiquark distributions according to the rule : 

u{-x) + d{-x) = -[u{x) + d{x)], (0<a;<l), (7) 

u{-x) - d{-x) = -[u{x) - d{x)], {0<x < 1), (8) 

An(-x) + Ad{-x) = Au{x) + Ad{x), {0 < x < 1), (9) 

Au{-x) - Ad{-x) = Au{x)-Ad{x), (0 < a; < 1). (10) 

Here, the long-dashed and dash-dotted curves respectively stand for the contribution of the 
discrete valence level (it is a deep bound state emerging from the positive energy continuum 
under the influence of hedgehog mean field) and that of the negative energy Dirac sea (the latter 
is sometimes called the vacuum polarization contribution), while their sums are shown by the 
solid curves. We emphasize that the separation into the valence and sea-quark contributions 
is highly model-dependent concept, and too much significance should not be given to it. Let 
us first look into u{x) + d{x) in Fig.l(a). This distribution, which emerges at the leading order 
of Nc, was first evaluated by Diakonov et al. |p. As emphasized by them, the proper account 
of the vacuum polarization contribution is essential here. In fact, the "valence-quark-only" 
approximation would lead to u{x) + d{x) < [|], thereby violating the positivity of the parton 
distribution. After including the vacuum polarization contribution, on the other hand, this 
fundamental requirement is satisfied finely. Next, Fig. 1(b) shows the isovector combination of 
the unpolarized distribution function u{x) — d{x). The vacuum polarization contribution is 
sizable also for this quantity. Among others, the significant positivity of u{x) — d{x) in the 
negative x region denotes that u{x) — d{x) < for physical value of x with < x < 1, thereby 
indicating the excess of d sea over the u sea in the proton. This result follows from the fact 
that the effect of pion cloud is automatically incorporated into the framework of the CQSM 
||23| . We shall later compare the above prediction of the CQSM directly with the existing high 
energy data. 

Next, we turn to the longitudinally polarized distributions. The isoscalar combination 
Au{x) + Ad{x) and the isovector one Au{x) — Ad{x) are respectively shown in Fig.l(c) and 
Fig. 1(d). The isoscalar distribution survives only at the first order of Q or equivalently at 
0{N^). This is because the leading-order term vanishes identically because of the hedgehog 
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symmetry. On the other hand, the leading contribution to the isovector distribution arises 
from the 0{fl^) or 0{N^) term. Although the main feature of the isovector longitudinally 
polarized distributions are controlled by this leading order contribution, the next-to-leading- 
order contribution is also important. In fact, without this 1/Nc correction, the first moment 
of Au{x) — Ad{x), i.e. the isovector axial charge of the nucleon g)^ would be largely un- 
derestimated. (Incidentally, the presence of this first order rotational correction to g)^ is a 
distinguishable feature of our effective quark theory as compared with a deeply connected 
soliton model of the nucleon, i.e. the Skyrme model |]2^,[^. This breakdown of the so- 
called Cheshire-cat principle can be understood as originating from the noncommutativity of 
the bosonization procedure and the collective quantization procedure of the rotational motion 
26|.) The numerical results also shows that the isoscalar and isovector distributions have 



rather different behaviors as functions of x. The first thing one may notice for the isoscalar 
distribution is the smallness of the sea quark contribution. (At least the smallness of the sea 
quark contribution to the first moment of Au{x) + Ad{x) has long been known to us 0.) The 
smallness of the sea quark contribution here does not necessarily mean negligible role of the 
vacuum polarization effect itself. Since our valence quark level is a deep bound state generated 
by the strong mean-field of hedgehog shape, it is legitimate to think that the valence level 
contribution above also contains some sort of vacuum polarization effect due to this strong 
mean field. In fact, the first moment of Au{x) + Ad{x), which can be interpreted as the quark 
spin content of the nucleon AS3 or the flavor- singlet axial charge at the energy scale of the 
present model, turns out to be fairly small as 

AS3 = J dx[Au{x) + Ad{x)] 

= dx[Au{x) + Au{x) + Ad{x) + Ad{x)] - 0.35. (11) 
Jo 

As advocated in 0, the symmetry breaking strong mean- field of hedgehog shape must be 
responsible for this smallness of the quark spin content. A quite interesting feature of the 
theoretical distribution Au{x) + Ad{x) is that it changes sign from positive to negative as 



X becomes smaller than about 0.1. As already pointed out in our previous paper [|I2|, what 
causes this sign change is the first order rotational correction, which arises from the proper 
account of the nonlocality (in time) of the quark bilinear operator entering into the theoretical 
definition of the quark distribution function. We shall later show that this sign change of 
the isoscalar distribution function in the small x region is just what is required by the recent 
experimental data for the deuteron structure function gf{x). 

The isovector combination Au{x) — Ad{x) illustrated in Fig. 1(d) shows totally different x 
dependence from the isoscalar one. The vacuum polarization contribution to this distribution 
function has a sizable magnitude of peak with positive sign around x ~ 0. A sizable nonzero 
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value of Au{x) — Ad{x) in the negative x domain has a far-reaching physical consequence. 
That is, it means the violation of SU(2) symmetry for the longitudinally polarized sea-quark 
(anti-quark) distributions. In consideration of the relation (|T0|), we can read from Fig. 1(d) 
that Au{x) — Ad{x) > for physical range of x between and 1. 

To sum up, within the theoretical framework of the CQSM, the isoscalar and isovector 
distribution functions turn out to have fairly different x dependence reflecting the fact that 
they have totally dissimilar dependence. Now we shall compare these predictions of the 
CQSM with the existing high energy data for the unpolarized as well as the longitudinally 
polarized structure functions. Before carrying out this comparison, we must first take account 
of the scale dependence of the distribution functions. We have done it by using the Fortran 
programs given by Saga group ||2^, which provide us with solutions of Dokshitzer-Gribov- 
Lipatov-Altarrelli-Parisi (DGLAP) evolution equation at the next-to- leading order (NLO). The 
question here is what value we should use for the initial energy scale of this evolution. Here 
we have tried to see the effect of variation of Qfnu ^ range centered around the model energy 
scale of Mj.y ~ (0.56 GeV)^ - 0.31 GeV^. We find that the final results are not significantly 
different in the range 0.25GeV^ < Qj^^ < 0.35GeV^. We shall use the value Qf^^ = 0.30 GeV^ 
throughout the following analyses. In solving the DGLAP equation, we set Nf = 3 and assume 
s{x) = s{x) = 0, g{x) = and As{x) = As{x) = 0, Ag{x) = at the initial energy scale. 



Figure 2: The predictions for -Ff (^) ~ ^-zi^) ^^'^ -F2"(^)/-^2 (^) = AGeV"^ are compared 

with the NMC data given at the corresponding energy scale |l|. 



Fig. 2 shows the theoretical predictions for -Ff (^) ~^ni^) -F"(a;)/-Ff (x) at = 4GeV^ 
in comparison with the NMC data. One sees that the NMC data for F2{x) — F^{x) are well 
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reproduced by the theory, except for a small discrepancy between the peak positions of the 
theoretical curve and the experimental data. Next, we compare our theoretical prediction for 
the ratio F^(x)/F2(x) with the corresponding NMC data. One can say that the agreement 
between the theory and the experiment is fairly good at least for x < 0.25. The theoretical 
prediction deviates from the data as x becomes larger than 0.3. This might be connected with 
the fact that our treatment of the soliton center-of-mass motion is essentially nonrelativistic. 
Note, however, that a reliable theoretical prediction of the ratio F^{x) / F^\x) near x ~ 1 is 
extremely difficult, because both of F"'{x) and F^[x) damps rapidly as x becomes larger. 



Figure 3: The theoretical predictions for the unpolarized antiquark distribution d{x) — u{x) 
at = 2.3 GeV and = 54 GeV in comparison with the HERMES H and the E866 data 
||T5| . Here the open squares represent the HERMES data given at = 2.3 GeV, while filled 
squares denote the E866 data given at = 54 GeV. 

Next, in Fig. 3, we compare the theoretical predictions for the anti-quark distribution d{x) — 
u{x) evaluated at two different values of with the HERMES |]14| and NuSea data |15[ given 



at the corresponding energy scales. (One certainly confirms that the scale dependence of the 
quantity d{x) — u{x) is rather small.) The general trend of the experimental data is well 
reproduced by the theory, although the theory has a tendency to slightly overestimate the 
magnitude of the fiavor asymmetry of the sea-quark distributions. 

Now we turn to the discussion of the longitudinally polarized distributions. Shown in 
Fig. 4 are the longitudinally polarized structure functions for the proton, the neutron and the 
deuteron, gi{x), g^i^x) and gf{x). The results for gi{x) and gi{x) were already reported in 
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T^, but we show them again for the purpose of comparison, since gf{x) is essentially given 
as their average. To be more precise, we evaluate gf{x, Q^) in Fig.4(c) by making use of the 
standardly-used formula 

gfi^, Q') = li 9l{x, Q') + Q')){1- 1.5un ), (12) 
with lod the probability that the deuteron is in a D-state. 



Figure 4: The theoretical predictions for the longitudinally polarized structure functions for the 
proton, the neutron and the deuteron at = 5 GeV^ in comparison with the corresponding 
experimental data. The filled circles in (a) and (b) respectively corresponds to the E143 |jl6| 
and the E154 data [|T7|, whereas the filled circles, the open circles and the open squares in (c) 
and (d) represent the E143, the El 55 fg and the SMC data [|1|]. 



As already argued in |T^, a prominent feature of the CQSM is the good reproduction 
of the neutron data, which are known to have large magnitudes with negative sign. We 
have interpreted this as a manifestation of the chiral symmetry, maximally incorporated into 
this effective quark model. Now a direct comparison with the deuteron data reveals another 
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(though related) interesting aspect of physics. An noticeable characteristic of the recent data 
for gf{x, Q^) is that it appears to show sign change as x approaches zero, although care must 
be paid to the fact that the precision of the experimental data in the small x region is still poor. 
Quite interestingly, the theoretical prediction of the CQSM closely follows the variation of this 
gf{x, Q"^) at least qualitatively. What is the origin of this behavior of the theoretical structure 
function? Undoubtedly, it can be traced back to the behavior of the isoscalar distribution 
function Au{x) + Ad{x) shown in Fig. 1(c). In fact, as emphasized by Windmolders |2^, the 



singlet distribution is strongly constrained by the measured value of gf. (This seems only 
natural, since the deuteron is an isoscalar target.) To see it more closely, we recall the general 
expression for gf given as 



. 1/1 



91 - q[ ^Cns ® Ags + C5 ® AS + 2iVjC, ® Ag j , (13) 
where 

Aqs = Au + Ad-2As, (14) 

AS = Au + Ad + As, (15) 

and Ciys,Cs, Cq are the non-singlet, singlet and gluon Wilson coefficients and the symbol ® 
represents convolution with respect to x. 

Assuming minor role of the s-quark degrees of freedom at the energy scale in question, the 
difference between Aqg and AS would be small. The gf is then determined by Am + Ad and 
Ag at this energy scale. Now, somewhat peculiar small-x behavior of the isoscalar longitudi- 
nally polarized distribution function predicted by the CQSM seems to get a phenomenological 
support of the latest deuteron data. In any case, the gf, especially in a smaller x region, 
appears to be extremely sensitive to the detailed dynamical content of the theories. To be able 
to reproduce it or not would then be a good touchstone of model selection. 

Now that we have confirmed that the CQSM can explain the principle experimental data 
for the unpolarized and the longitudinally polarized structure functions of the nucleon as well 
as the deuteron, we come back to an interesting prediction of the CQSM given by Fig. 1(d), i.e. 
the possible violation of the light fiavor sea-quark asymmetry for the longitudinally polarized 
distributions. Since the sea-quark distributions are not well constrained by the inclusive data 



alone, very little is known as to them. In a recent report [|29[, Morii and Yamanishi proposed 
new formulas for extracting a difference, Ac? — Au, from data of the polarized semi-inclusive 
processes / + A^^ I' + H + X and estimated the value of it from the available data of the 
SMC and Hermes groups. The results of their analyses are shown in Fig. 5 together with the 
prediction of the CQSM at the corresponding energy scale, = 4GeV^. Unfortunately, the 
uncertainties of the fit are so large that nothing definite can be said at the present status. 

10 



Figure 5: The theoretical prediction for the longitudinally polarized antiquark distribution 
Ad{x) — Au{x) at = 4GeV^ is compared with the recent semi-phenomenological fit by 
Morii and Yamanishi. See figure caption of |^ for the detailed meaning of three different 
marks in their results. 

Note, however, that Morii and Yamanishi parameterized Ad{x) — Au{x) as 

A - Au{x) = C ( d{x) - u{x) ) (16) 

and determined the value of C and a from the x^-fit of the results presented in the above 
figure. The results were C = —1.00 and a = 0.18, thereby indicating an asymmetry of d and 
u. The negative value of C denotes that 

Ad{x) - Au{x) < 0, (17) 

since we already know that d{x) —u{x) > 0. This is qualitatively consistent with the prediction 
of the CQSM. In fact, the values of C and a determined from our theoretical distributions at 
= 4 GeV^ turn out to be 

C--2.0, a -0.12, (18) 

in qualitative consistency with the result of their analysis. 

Despite encouraging phenomenological support to the present approach, an important ques- 
tion remains to be answered. The question concerns the role of gluons at the relatively low 
energy scale of Q2 = o.3 ~ 0.5 GeV^ which may be taken starting energy of the DGLAP 
evolution equation. We have simply assumed zero for the input gluon densities at the initial 
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energy scale. This assumption may not be necessarily justified, especially for the unpolarized 
case. In fact, there is some phenomenological indication that the gluon saturates nearly 30 % 
of the nucleon momentum sum rule even at such low energy scale (Note, however, that 
the isovector-type quantities like -Ffl^) ~ F2{x) and d{x) — u{x) investigated in the present 
paper would be insensitive to these input gluon densities.) On the other hand, very little is 
known for the polarized gluon densities. A relatively good agreement between the experimen- 
tal structure functions and our theoretical predictions, which are obtained by assuming zero 
gluon polarization at Qfnu = 0.3 GeV^, indicates that the role of gluons may be less important 
than the case of the unpolarized distributions. (This never denies the importance of gluon 
polarization at high energy scale, since the gluons are known to acquire polarization rapidly 
through the processes of perturbative evolution.) Undoubtedly, it is a very important the- 
oretical question to reliably predict the gluon densities to be used as input distributions of 
the renormalization-group evolution. The importance of the gluon densities at this matching 
energy scale of = 0.3 ~ 0.5 GeV^ was also emphasized by by Lampe and Reya in a recent 



review According to these authors, if low energy models cannot provide necessary input 
gluon densities at this energy, they would only refer to some nonperturbative input quark scale 
which cannot be reached by perturbative evolution. 

In summary, the CQSM allows us for the calculation of the full x-dependence of the nucleon 
structure functions without introducing any new free parameters. In particular, what makes 
this approach quite promising is the fact that it predicts, besides the valence-like distributions, 
the sea-quark-like densities concentrated in the small x domain, if it is treated in a theoretically 
consistent manner (i.e., not in the "valence-quark-only" approximation). This unique feature of 
the model is expected to provide us with a good starting point for theoretically understanding 
unpolarized data as well as yet-to-be-obtained polarized data in the small x region. Making 
full use of this advantage of the CQSM, we have made a prediction for the flavor asymmetry of 
the longitudinally polarized antiquark distributions in the proton. The model has been shown 
to predict Ad{x) — Au{x) < 0, while at the same time d{x) —u{x) > 0, i.e., it definitely predicts 
light flavor sea-quark asymmetry not only for the spin-averaged distributions but also for the 
longitudinally polarized distributions. We hope that this interesting prediction of the CQSM 
will be tested through accumulation and analyses of more precise experimental data on the 
polarized semi-inclusive processes in the near future. 
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